Anomalous spin precession and spin Hall effect in semiconductor quantum wells 
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Spin-orbit (SO) interactions give a spin-dependent correction r so to the position operator, referred 
to as the anomalous position operator. We study the contributions of r so to the spin-Hall effect 
(SHE) in quasi two-dimensional (2D) semiconductor quantum wells with strong band structure SO 
interactions that cause spin precession. The skew scattering and side-jump scattering terms in the 
SHE vanish, but we identify two additional terms in the SHE, due to r so , which have not been 
considered in the literature so far. One term reflects the modification of the spin precession due 
to the action of the external electric field (the field drives the current in the quantum well), which 
produces, via r so , an effective magnetic field perpendicular to the plane of the quantum well. The 
other term reflects a similar modification of the spin precession due to the action of the electric field 
created by random impurities, and appears in a careful formulation of the Born approximation. We 
refer to these two effects collectively as anomalous spin precession and we note that they contribute 
to the SHE to the first order in the SO coupling constant even though they formally appear to 
be of second order. In electron systems with weak momentum scattering, the contribution of the 
anomalous spin precession due to the external electric field equals 1/2 the usual side-jump SHE, 
while the additional impurity-dependent contribution depends on the form of the band structure 
SO coupling. For band structure SO linear in wave vector the two additional contributions cancel. 
For band structure SO cubic in wave vector only the contribution due to external electric field is 
present, and can be detected through its density dependence. In 2D hole systems both anomalous 
spin precession contributions vanish identically. 



I. INTRODUCTION 

In systems with strong spin-orbit (SO) interactions 
an electric field generates a transverse spin-current 1-12 : 
this phenomenon is referred to as spin-Hall effect (SHE). 
For the past ten years, the SHE has been a source 
of new ideas for magneto-electronic devices 13 aimed at 
integrating semiconductor and magnetic technologies, 
facilitating efficient information processing and quan- 
tum computing architectures. 14-17 These visions have 
stimulated a large volume of experimental and the- 
oretical work. 18-33 Experimentally, the SHE was ini- 
tially studied in semiconductors, 34-38 but has since ex- 
panded to novel materials such as HgTe-based quantum 
wells, 39 topological insulators and graphene, and d-band 
metals. 40-42 It is often simpler to measure the inverse 
spin-Hall effect, 43 where a spin current generates a trans- 
verse charge current, which is detected by conventional 
means. The inverse SHE has been observed in Al, 44 Pt 
wires at room temperature, 45 hybrid FePt/Au devices, 46 
Au films with Pt impurities, 47 permalloy/normal metal 
bilayers, 48 GaAs multiple quantum wells, 49 and Cu with 
Ir impurities. 50 For a review of recent experimental work 
on the SHE in Pt see Ref. 51. Observation of the in- 
verse SHE has recently been reported even in a weakly 
spin-orbit coupled material such as Si. 52 

SO coupling may be present in the band structure and 
in the impurity potentials. Band-structure SO interac- 



tions become important in structures lacking a center of 
inversion when SO interactions lift spin degeneracy. 77 If 
the underlying crystal lattice lacks a center of inversion 
the material is said to possess bulk inversion asymme- 
try (BIA). In low-dimcnsionsional systems the confine- 
ment potential can be made asymmetric, in which case 
one speaks of structure inversion asymmetry (SIA). In 
this paper we consider exclusively quasi two-dimensional 
semiconductor systems that lack a center of inversion due 
to BIA 78 and/or SIA giving rise to Rashba SO coupling. 
In these systems the band structure SO interaction is 
represented by a Hamiltonian H = (h/2) cr ■ $7*. describ- 
ing the interaction of the spin with an effective wave 
vector-dependent magnetic field fife . This can be fl^ 1A or 
r?| IA . The spin precesses about this field with frequency 
ilfe = | fife |- Different physical regimes are distinguished 
by the value of the product of f2fc with the momentum 
relaxation time r p . In the ballistic regime f2fcT p — > oo. 
The weak momentum scattering regime is characterized 
by QkTp 3> 1, while in the strong momentum scattering 
regime f2feT p <C 1. 

SO interactions arise, quite generally, from a spin- 
dependent correction r so to the position operator, 53,54 
whose general form is 

f so = Ao-xwfe, (1) 

where A is a material-specific parameter, cr is the vector 
of Pauli spin matrices, and u)k takes different forms for 
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different systems, as weil as for electrons and holes in the 
same system. 79 We note that fl^ and us^ are not inde- 
pendent of each other. The uJk entering the corrected 
position operator is inherently related to J~2| IA charac- 
terizing the Rashba SO coupling through a term of the 
form $~2| IA = \u>k x W. The presence of r so results 
in corrections to the interaction between charge carriers 
and electric fields, which include impurities and external 
electric fields. Thus, in addition to the band structure 
SO interaction, one must take into account SO interac- 
tions arising from the external electric field and from the 
electron-impurity potential. The interplay between these 
interactions in the SHE is quite a complicated subject. 
It has received a lot of attention in recent years, yet, as 
we will see, it is not yet completely understood. 

Perhaps the most intuitive mechanism of SHE is the 
one known as skew scattering, i.e., the asymmetric scat- 
tering of up and down spins by impurities. 55-57 Next, 
we have the so-called side-jump scattering 53 ' 58 ~ ea which 
consists of two equal terms, one reflecting the correction 
to the band energy due to the spin-dependent interaction 
with the electric field, the other reflecting the renormal- 
ization of the carrier trajectory during collisions. Dia- 
grammatic formulations naturally recover the two side- 
jump scattering terms through the vertex renormaliza- 
tion of spin and charge currents, as Ref. 56 demon- 
strated. An analytical derivation of the side jump from 
the Kubo formula was presented in Ref. 61. Further- 
more, Ref. 10 identified skew-scattering and side-jump 
scattering within a drift-diffusion approach. More re- 
cently, side jump scattering was derived starting from the 
quantum Liouville equation for the single-particle spin 
density-matrix. 62 

The analysis of the SHE becomes considerably more 
complicated when both band-structure and impurity- 
potential induced SO interactions are present. This prob- 
lem was first addressed by Tse and Das Sarma, 63 who 
employed the diagrammatic Kubo formula and consid- 
ered band structure SO coupling of the linear Rashba 
form. They found that the skew scattering contribution 
to the SHE vanished for arbitrarily small value of the 
band structure SO coupling, while a term equal to half 
the usual side-jump scattering SHE survived. 80 This is 
in contrast to the result obtained in Ref. 64 that both 
the side-jump and the skew scattering contributions van- 
ish for arbitrarily small value of the band structure SO 
coupling, as long as the impurity-induced (Elliot- Yafct) 
spin relaxation is neglected. These two results are recon- 
ciled by taking into account the SO contribution to the 
electron-impurity self-energy diagram, 65 which recovers 
the vanishing of the side-jump and skew scattering con- 
tributions found in Ref. 64. 

The principal question identified in Ref. 63 was the 
paradox of the non-analyticity of the spin Hall conduc- 
tivity, which appears to change discontinuously as soon as 
the band structure SO coupling is turned on. This para- 
dox was finally solved in Ref. 66 by the introduction of 
an impurity-induced (Elliott- Yafet) spin relaxation rate 



1/teY; which led to a spin Hall conductivity of the form 

, [ <J yx}ss + [0'yx}s3 , . 

a yx = -i T ~ JZ y Z > 

1 + tey/tdp 

where tdp is the Dyakonov-Perel relaxation time asso- 
ciated with the band structure SO coupling and given 
by = (0 2 )t p , where r p is the momentum relax- 
ation time and the angular bracket denote an average 
over the momentum distribution. The above formula 
exhibits a smooth crossover between the sum of skew- 
scattering (ss) and side-jump (sj) contributions, when 
the band structure spin precession VL^ is neglected, and 
zero when f2fcT p 3> 1, i.e. when the band-structure SO 
interaction is much stronger than the electron-impurity 
interaction (see also Ref. 65). 

However, this is not the end of the story. The work 
described above was limited to band-structure SO cou- 
plings that are linear in wave vector k. The aim of this 
work is to provide a consistent framework for treating 
band structure and impurity SO effects in quasi two- 
dimensional quantum wells for any form of the band 
structure SO interaction in the weak momentum scat- 
tering regime f2fcT p 3> 1. To this end, we construct a 
rigorous theory of the interplay of spin precession due 
to band structure SO coupling and SO coupling due to 
impurities. We start from the quantum Liouville equa- 
tion and derive a kinetic equation for the spin density 
matrix, which captures the effects of band-structure spin 
precession and r so on an equal footing. We focus from 
the very beginning on the weak momentum scattering 
regime VLf.r p 3> 1. Under this assumption, we do not 
have to worry about the finite Elliot- Yafet scattering rate 
that appears in Eq. (2): we are in the regime tey 3> tdp- 
But, while Eq. (2) predicts, in this limit, a vanishing spin 
Hall conductivity for linear-in-fe band-structure SO inter- 
action, we will show that a finite spin Hall conductivity 
can survive for different forms of that SO interaction. 

More precisely, we find that, in the weak momentum 
scattering regime, skew scattering and side jump scat- 
tering still give zero SHE. At the same time, we iden- 
tify two additional contributions to the SHE stemming 
from r so . These contributions have been overlooked in 
the literature thus far. One contribution arises from the 
impurity potential, and is found in the Born approxi- 
mation when scattering terms of second order in SO are 
taken into account. This contribution can be viewed as a 
modification of the band structure precession frequency 
due to electron-impurity interaction. The second con- 
tribution is scattering- independent. Its origin lies in the 
spin-dependent interaction with the external electric field 
brought about by f so . This looks like each carrier inter- 
acting with an effective magnetic field. The carrier spin 
precesses in this effective magnetic field in such a way 
that an out-of-plane spin component is generated, which 
contributes to the SHE. We refer to these two effects 
collectively as anomalous spin precession. The impurity- 
induced anomalous spin precession term gives an out-of 
plane component of the effective magnetic field. This 
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is precisely what distinguishes the anomalous spin pre- 
cession from the usual side-jump scattering term, which 
vanishes in the presence of spin-precession. Remarkably, 
these effects contribute to the SHE to the first order in 
the SO coupling constant even though they formally ap- 
pear to be of second order. The external electric field part 
of the anomalous spin precession appears to be universal 
in electronic systems in the clean limit. 

In electron systems with band structure SO linear in 
k the sum of the two anomalous spin precession terms 
vanishes. In hole systems both additional terms arc zero 
independently. Nevertheless, the electric field part of the 
anomalous spin precession term in the SHE in general 
survives, and we demonstrate its existence explicitly in 
2D electron systems with band structure SO described 
by the cubic Dresselhaus model. This model applies to a 
wide quantum well at high density n e , in which we find 
the total SHE conductivity in the clean limit to be 



The term oc A is detectable since it is linear in the electron 
number density, while the band-structure SO contribu- 
tion in the weak momentum scattering regime is density- 
independent. Our results are relevant to experiment and 
help to distinguish different contributions to the SHE. 

We do not discuss explicitly contributions to the SHE 
purely from band structure SO, which are well known. 6 
The focus of this work is on the contributions to the 
SHE due to r so , and the central result is that, aside from 
the well-known skew scattering and side-jump scatter- 
ing terms, two additional contributions - the anomalous 
spin precession terms - are present when band structure 
SO is nonzero. To our knowledge, this is the first work 
that proves that r so can give rise to a spin-Hall current 
through a mechanism unrelated to scattering. We work 
up to third order in the impurity potential, and in or- 
der to recover all contributions consider terms of second 
order in the SO coupling. Our results are valid in the 
weak momentum scattering limit, yet in the Appendix 
we prove rigorously that a non-analyticity in the strong 
momentum scattering limit is cured by introducing the 
Elliott- Yafet spin relaxation time tey, as was done in 
Ref. 65. 

The outline of this paper is as follows. In Sec. II we 
present the band Hamiltonian and in Sec. Ill we dis- 
cuss the effective position operator. In Sec. IV we derive 
the general form of the kinetic equation starting with the 
quantum Liouville equation, and discuss the various scat- 
tering terms. In Sec. V we discuss the non-equilibrium 
correction to the density matrix, demonstrating that a 
new, scattering-independent driving term due to r so is 
present. The general solution to the kinetic equation is 
presented in Sec. VI, demonstrating that the skew scat- 
tering and side-jump scattering terms give zero contribu- 
tions to the SHE. All SHE contributions due to r so are 
listed for commonly employed models of SO coupling. 
An explanation of anomalous spin precession is given in 



Sec. VII, which is followed by a summary and conclu- 
sions. 



II. BAND HAMILTONIAN 

In the crystal-momentum representation, the band 
Hamiltonian Hq in the effective mass approximation has 
the general form 

Hok = -ffkin + H so = i?km + - er ■ Ofc, (4) 

for an arbitrary SO interaction. The kinetic energy term 

#kin = £ok 1 = f-^r 1 , where 1 is the identity matrix 
in spin space and m* the carrier effective mass. The 
spin-dependent term in the Hamiltonian H so is treated 
as a perturbation with respect to the kinetic energy term. 
The eigen-energies are written as Ek± = £ok ± {Klk/2). 

For 2D spin-1/2 electron systems with SIA, the band 
structure contains the linear Rashba Hamiltonian 

Hri = oti (a x k y - <T y k x ) = aii{k_a + - fc + er_), (5) 

where k± = k x ±ik y and a± = (a x ±ia y )/2. In these sys- 
tems BIA has two contributions, the linear Dresselhaus 
term 

Hdi = Pi(cr x k x - dyky) = /3i(fc+(7 + + /s-cr-), (6) 
and the cubic Dresselhaus term 
#153 = MvxKkl - Vykykl) 

= h[k-(k\ - k 2 _)a+ + k+{kl - k\)a-]. (7) 

For 2D heavy-hole systems SO coupling due to SIA is 
dominated by the cubic Rashba Hamiltonian, 

Hr3 = a 3 [k y (k 2 - 3kl)a x + k x {k 2 x - k\)a y \ 

= otii(k\o ■_ - kta +) . (8) 

BIA in 2D heavy-hole systems contains the fc-linear term 
H D v = li{°xk x + Oyk y ) = 7i(fc + CT_ + k-a + ), (9) 
and the cubic Dresselhaus term 

Hd3' = lz{kl + ky)(a x k x + a y k y ) 

= l3 (k 2 + k_a- + k 2 _k + a+). (10) 

For the terms cubic in k, we restricted ourselves to the 
dominant contributions due to SIA and BIA. 

III. EFFECTIVE POSITION OPERATOR 

The SO interaction appears when transforming from 
the Dirac to the Pauli equation by means of the Foldy- 
Wouthuysen transformation. 54 Under this transforma- 
tion, the position operator in spin-1/2 systems becomes 

f p hys = r + r so , (11) 
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where the SO part r so is expressed in terms of the vector 
a of Pauli spin matrices. We refer to r so as the anoma- 
lous position operator. 

The general form for the anomalous position operator, 
valid for both 2D electron and 2D hole systems, is 



Xcr x o> fc , 



(12) 



where A and u> k are different for electrons and holes. For 
2D electrons u) k = k, and 



r so = Ai <t x k, 



(13) 



assuming Xik F -C 1. For 2D hole systems the correction 
to the position operator has the form 



r so = A 3 er x u k3 , 



(14) 



where uj k3 = k 3 (cos 36*, sin 36*, 0), assuming X 3 k F <C 1. 

Consider a general scalar potential V(f). Under 
the Foldy-Wouthuysen transformation it transforms to 
V(f p hys), which, to first order in r so , takes the form 

V(r phys ) = V(r) + 1 [W(r) • r so + f so • W(f )]. (15) 

Therefore, as a result of this transformation, both the 
potential due to an applied electric field and the impurity 
scattering potential acquire spin-dependent terms. 

Let U(r) denote the scattering potential, which repre- 
sents elastic scattering off charged impurities and static 
defects (but not phonons or electrons) 



U(r)=J2u(r-R r ), 



(16) 



where Ri indexes the random locations of the impurities 
and the scattering potential due to a single impurity is 
denoted by U(r). In Fourier space, the matrix elements 
of U(r) are 



(17) 



and the potential due to a single impurity is written as 

Ukk> =U kk ,t + V kw , (18) 

where lA kk i represents the matrix element of the poten- 
tial due to a single impurity between plane waves, while 
Vfcfc' is the spin-dependent part arising from r so . Both 
have units of energy x volume. The strength of the dis- 
order potential is characterized by the impurity density 
rii. The matrix elements of the spin-dependent part of 
the impurity potential in reciprocal space are 

iX 

V kk > = -— a ■ (w fe x k' -u k > x k)U kk >. (19) 

In 2D the spin dependent term in V kk i points out of the 
plane for both electron and hole systems. 



Interaction with a static, uniform external electric field 
E is contained in 

H E kk' = (eE ■ r) kk 't + e (E ■ r so ) kk S kk ' 

d 1 (20) 

= ieE ■ — S(k - k') t + - a ■ A k 5 kk , . 

with 1 the identity matrix in spin space, and arises 
from the anomalous position operator. 67 From Eq. (12), 



Afc = 2eAwfc x E. 



(21) 



It follows from the preceding discussion that has dif- 
ferent forms in electron and hole systems. 

The anomalous position operator accounts for impu- 
rity SO coupling and for band structure SO coupling due 
to SIA. To see the latter, consider the SO coupling due 
to the full potential Vtot acting on the system. In a 2D 
system we can divide Vtot = V ex t + Vqw + U , where V ex t 
is the applied electric field, Vqw the z-direction confine- 
ment, and U the impurity potential introduced above. 
The total potential Vtot gives rise to SO coupling, which 
in reciprocal space is contained in 

H so ,k = X n a- ■ k x V(V ex t + U) + X n er -kx z {^—^ 

Z (22) 

In the second term we can incorporate the average 
(9Vqw/9z) over the quantum well into an effective SO 
constant a, giving the Rashba SO coupling. 81 This clar- 
ifies the relationship between a and A and shows that, 
knowing the form of the Rashba Hamiltonian in a certain 
system, one can deduce the form of r so in that system. 

The full Hamiltonian is H k ot = H ok + HE kk > + U kk i . 
The spin current operator jj corresponding to spin com- 
ponent i flowing in the direction j is 



H 2 kj 
2m 



(23) 



In addition to the contribution from the band Hamil- 
tonian, the velocity operator has two additional terms, 
discussed in detail in Ref. 62. The first stems from the 
spin-dependent interaction with the external electric field 
HEk\, while the second arises from the spin-dependent 
term V kk i in the impurity potential. These two cancel, 
as they represent the net force acting on the system. 62 
They will not be explicitly considered in what follows. 



IV. KINETIC EQUATION 

The formalism presented here parallels that originally 
formulated in Rcfs. 8,9. The Liouville equation for the 
density operator p is projected onto the basis {|fe)}, with 
Pkk' = fk Skk' +9kk>, where g k k> is off-diagonal in fc, and 
all quantities are matrices in spin space. These satisfy 

dfk 



(24a) 



dgkk' 
dt 



[Ho,9] 



kk' 



[UJ}kk> -r-[U,g} kk ,,(24b) 
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We focus on variations which are slow on the scale of 
the momentum relaxation time, and solve for g kk > as an 
expansion in the impurity potential, which can be per- 
formed to any desired order. Very generally f k satisfies 

^ + ^[H o J] kk + J(f k ) = 0. (25) 

The total scattering term J(f k ) = JBom(fk) + J ss (fk), 
where in the first Born approximation 

■W/fe) = ^(J^^ ft, e- i6ot '/ n [U, f] e^'V^A } 

(26) 



and (...) represents averaging over impurity configura- 
tions. In the second Born approximation we obtain the 
additional skew scattering term 



/ />oo />oo V 

J ss (fk) = -T E { dt' rfi , q^,e-' i ^ t '/ fi [L>,e- ldot ''/ n [L>,/]e^ ot '' /n ]e^ ot '/ ?l ] ) . (27) 

* \ J0 JO I kk 

We expand JBorn(/fe) in ftk and A. We retain the leading term plus terms to first order in fife, first order in A, and 
the second-order term in fifcA. Thus JBom(/fe) can be written as a perturbation expansion in tt k and A in the form 

J B orn(/fe) = Mfk) + Mfk) + Mfk) + Jnx(fk)- (28) 

The leading term in JBom(/fc) is the scalar Jo(/fe), which is the customary Born-approximation scattering term 
appearing in the Boltzmann equation. It is found by taking Eq. (26) and considering only the scalar parts of Hq (i.e. 
-ffkin) and U (i.e. Ukk'), and in 2D takes the form 

Mfk) = J J \Ukk>\ 2 {fk fk>). (29) 

Next, we have the term in JBom(/fc) to first order in fife (i.e. due to band-structure SO coupling), which is found by 
considering the spin-dependent part of Hq and the scalar part of U . It gives rise to a well-known scattering term, 
referred to here as Jo(/fc). 8 ' 9 We only require its action on the scalar part of the density matrix, rife, given by 

7r f d 2 k' d 
Jq. (rik) = t / 77, — |£4fc' I 2 ("fe ~nk')cr- (fl k - Si k >) it— 5{e ok - e fc')- ( 30 ) 
h J (2ttY de 

This term is relevant only in determining the band-structure SO contribution to the spin current, which has been 
studied previously, and is not pertinent to the discussion presented in this work and will not be given. Following 
on, in the side-jump scattering term J S j(f k ) we take the scalar part of Hq and the spin-dependent part of U. The 
electric field E is also finite in this term: without it J S j{f k ) would vanish. 82 Because E is nonzero, J S j(f k ) acts on 
the equilibrium density matrix fo k . It has two parts, which have been determined in Ref. 62. We use the notation of 
Ref. 62. We write J s j (n k ) = {n k ) + (n k ). The first part of the side-jump scattering term, referred to as J S j (n k ), 
arises from the change in the band energy due to the spin-dependent energy of interaction with E 

2im f d 2 k' 1 d 

(nfe) = — r- 1 -7— 2 \U kk ,\ 2 (n k -n k ,)-cr-(A k -A k ,)- S(s ok -s ok ,). (31) 

J h J (2tt) 2 2 ds ok 

The second part, J^. (n k ), reflects the spin-dependent change in the carrier position during collisions 



% i iriiTreE f tfk' ( ' dV k > k , dVy k 

sjv h J (27r) 2 ""™ V 9k' ' dk J v ~" '' de ok > 



Js] (nk) = ~z / , , U kk > H — — (rife - n k >) ^ S(e ok - e ok >) + h.c. (32) 



Both parts of the side jump scattering term are oc <j z . 
The scattering term Jn\{n k ) reads 



" 7TTL' /" d k! d 

Jn\(n k ) = / , n NW [a- ■ tt k > , V kk <] U kk > {n k - ny) — <5(eofc - eofc')- (33) 

h J (2ir) a de ok 



r 



The physical meaning of this term is as follows. During spin direction, given by fife, which represents the band- 
a scattering process, an incoming spin has a well-defined 
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structure SO coupling at wave vector k. Because the 
scattering potential is also spin dependent, the incom- 
ing spin is rotated during scattering by an amount that 
is proportional to Vkk' , the impurity SO coupling. This 
scattering term therefore represents spin rotations dur- 
ing collisions induced by the impurity SO coupling, the 
rotation being evident from its commutator structure. 

Even though we are doing perturbation theory to first 
order in the SO interaction terms A and ftk, spin preces- 
sion makes it necessary to include driving terms to order 
Af2fc, since these terms also yield contributions to the spin 



current oc A only, i.e. to first order in the impurity SO 
coupling. The necessity of including terms oc Xflk will 
become apparent when we discuss explicitly the solution 
for Seic introduced below, during which it will emerge 
that spin precession introduces a factor of 1/ ilfc . 

Beyond the first Born approximation we retain the 
leading term J ss (fk), in which A is finite but the elec- 
tric field E — 0, which is customarily responsible for 
skew scattering. 67 To first order in A, the real part of 
this term reduces to 



Jss{nk) = - 



37r 2 n,A f d 2 k 



d 2 k h 



Ukk'U k 'k"U k "k cr ■ (wfc x k! - u> k > x k) (rife/ - n k ")6(e k ~ £ofe")<H £ ofc - eofe')- 

(34) 



In 2D systems, in which both k and oj k are in the xy-plane, the skew scattering term is oc a z 



V. NON-EQUILIBRIUM DENSITY MATRIX 



In a constant uniform electric field E the density matrix is = /ofc + fEk- The equilibrium density matrix is given 
by 

/ofc = \ [/FD(£fe+) + /FD(£fe-)] + \ [/FD(£fc+) - /FD^fe-)] ° ' &k, (35) 

with /fd the Fermi-Dirac distribution function, while fEk is due to E. To first order in E the correction fEk satisfies 



g^ fc + ^ [H k , fEk] + J (fEk) 



eE dfpk 
h ' dk 



(36) 



r 



The term (eE/h) ■ (dfok/dk) corresponds to the usual 
streaming term in the Boltzmann equation. The second 
term on the RHS of Eq. (36) appears due to the anoma- 
lous position operator and is oc A. 

We write fk = rifcl + S'fc, where Sk is a 2 x 2 Hcr- 
mitian matrix, and correspondingly fEk = n-Ek^ + SEk 
and /ofe = n-ofcl + Sofc. The expectation values of the 
spin current operator is found from SEk- The term 
(eE/h) ■ (dfok/dk) may be decomposed into a scalar part 
(eE/h) ■ (dn k/dk) and a spin-dependent part (eE/h) ■ 
(dSok/dk). The spin-dependent part has been studied 
previously, 8,9 and is responsible for current-induced spin 
polarizations and spin currents arising from the band- 
structure SO coupling. It will not be discussed in this 
work. 



The non-equilibrium correction to the scalar part of 
the density matrix, nEk, is determined from 

dn E k , $ , \ eE) dn Qk 

-dT + Jo{nEk) = — -^k- (37) 

The solution to this equation is well known, and reads 
nEk = (&Et p /K) ■ (dnak/dk), with t p the momentum 
relaxation time. Once this solution is found, the spin- 
dependent scattering terms J ss , J S j and Jq\ act on nEk 
and produce additional effective driving terms for SEk- 
(The method used is the same as in Ref. 68.) 

We seek the solution for SEk to first order in A which 
we denote by SEkx- Specifically, including the contribu- 
tion due to Afc from Eq. (36), it is found from 



— ^ [Hk,SEk\\ + Jo(SEk\) = -J ss (nEk) - J S j(nEk) - Jn\(nEk) - ^ [Hsk\, Sok]- (38) 

I 

We specialize to short-range impurities henceforth, with- in Fourier space is written as 

out loss of generality. The potential of a single impurity _ ^ ^ (39a) 

iXU 

Vkk' = — a- ■ (uj k x k' - uj k ' x fe), (39b) 
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where the Fourier transform Ukk' has become the con- 
stant U. We write Jo(fk) = (/fc ~~ fk)/i~, with the over- 
line denoting an angular average over the directions of k, 
which in 2D indicates an average over the polar angle 9, 



X = 



dO 
2^ 



X, 



(40) 



and the momentum relaxation time tv, 



given by 



nim*U 2 



(41) 



We discuss the driving terms in more detail. Firstly. 

in.Xm* 2 ^ 



-J S s{n E k) 



a ■ (ojk x k' - u: k ' x k) n E k', 

(42) 



where the overlinc denotes averaging over 6' and the in- 
tegration over k' forces k' = k. We have established that 
this term is oc a z , and inspection of Eq. (42) reveals that 
this term is an odd function of k. 

The anomalous interaction with E gives rise to two 
driving terms. The first arises from the side-jump scat- 



tering term, which was determined in Rcf. 62. For both 
electrons and holes this takes the form 



Jsj( n Ek) 



a ■ A k S(e k - e F ) 



(43) 



This term is also odd in k. An additional driving term 
comes from the commutator of | a ■ with the den- 
sity matrix. Given that is already first-order in E 
we require only the equilibrium density matrix /ofc. We 
expand f k = fFD(e k )t + (h/2)cr ■ fi 



!FD{£k)lL "I" \nj L) <7 • \ Lk — 

the first term is a scalar, and at temperature T 
can write 



a/FD(£fc) where 
= we 



S3 



L'+Ehi foh] = s ( £ ok - £f) cr ■ rife x Afe. (44) 



Notice that this term is zero in the absence of spin preces- 
sion, when f k is a scalar and the commutator vanishes. 

The remaining driving term is —JnxijiEk)- For 2D 
electron systems, 



J 



? nr\ ni \U 2 f d 2 k' , d 

JnxynEk) = r / , n , 9 [cr ■ ilk' ,er ■ k x k \ [risk ~ n E k>) t\ o(s k - £ok')- 

n J [Zir) z OEQk 



(45) 



For 2D hole systems, 



- iXTrnAUl f d 2 k . d . ,.„. 

Jn\{n E k) = — / ,„ . 2 [cr flk>,cr ■ (w fc3 x k - u) k 'z x k)\ (n Ek - n Ek >) -= d(e fc - eofe'j- ( 46 ) 

2H J {2ny de k 



r 



VI. SOLUTION OF THE KINETIC EQUATION 



We summarize first the general solution to the kinetic 
equation for short range impurities and weak momentum 
scattering. We denote the driving terms generically by 
T^Ekx in this section. Let the component i of the spin 
operator be denoted by Sj = (h/2) <Tj. The spin density 
is Tr p§i = Tr psi , where the overbar denotes an angu- 
lar average as above, thus p is the isotropic part of the 
density matrix. Similarly, the spin current operator j'j 
has been defined in Eq. (23). Because it is odd in k 
its expectation value yields Tr pjj = Tr (p — p)jj . Con- 
sequently, the isotropic part of the spin density matrix 
determines the spin density, while the anisotropic part 
of the density matrix determines the spin current. It is 
therefore convenient to divide the spin density matrix 
into S Ek \ = S Ek \ +T Ek \, the isotropic part being S Ek \ 
(which gives the spin density) and the anisotropic part 
T E kX (which gives the spin current). From the quantum 
Liouvillc equation, we obtain a set of coupled equations 
for S E k\ and T E k\ for short-range impurities, which are 



solved rigorously in Appendix A. Here we just quote the 
solution. Letting T> E k\ = — cr d E k\, we find for fi^r ^> 1 



rp i (n h 
TEkX - 2 a ■ \JTk 



[d 



EkX 



- A 1 (d Ek x - Ad Ek \)}, 



(47) 

where the (dimensionlcss) matrix A is given by Aij — 
(Sij — fijOj), and T E k\ as found in Eq. (47) gives the spin 
current in the weak momentum scattering limit. Finally, 
we take the electric field E || x, the spin-Hall conduc- 
tivity is defined by — <j z yx E x , and we abbreviate the 
spin-Hall conductivity due to T E k\ simply by <j\. 



The appearance of the Qk in the denominator of Eq. 
(47) is a crucial feature of this solution. It demonstrates 
the need to retain scattering terms oc Afifc that are for- 
mally of second order in the SO coupling. 
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Table I: r so contributions to the SHE in units of n e eA for 
Qt S> 1. Here e~ stands for electrons and h + for holes, while 
Band SO abbreviates Band — structure SO. 



System 


Band SO 


A fe 


prcc 

a x 


< 


e~ 


Rl 


fc 


1/2 


-1/2 


e~ 


Dl 


fc 


1/2 


-1/2 


e~ 


D3 


k 


1/2 


-1/160 


h + 


m 










h + 


Dl' 


k" 








h+ 


D3' 


fc 3 









A. Skew scattering and side-jump scattering 



where n e is the electron density. In the weak momentum 
scattering limit this term is also independent of t. In 
2D electron systems it recovers the nonzero contribution 
to the SHE originally found by Tse and Das Sarma 63 
and subsequently by Raimondi and Schwab. 66 In 2D hole 
systems it is easy to check that cr^ rcc = 0. 

The origin of this contribution to the SHE will be elu- 
cidated in Sec. VII, but one remark is in order here. The 
spin-Hall conductivity cr^ rcc found in Eq. (49) has the op- 
posite sign to that found in Refs. 62,63,66 for the same 
orientation of the electric field. One should therefore not 
think of c^ lcc as a surviving side-jump term, but a quali- 
tatively new term due to r so altogether, which we identify 
with a spin precession mechanism with no counterpart in 
systems without band structure SO coupling. 



We recall that, as shown in Eqs. (42) and (43), both 
J ss (risk) and J s j (n F k) are odd in fc. Therefore the driv- 
ing terms due to J ss (riEk) and J S j (nEk) yield correc- 
tions to SEk\ that are even in fc. Since the spin current 
operator jj is odd in k, simple power counting in Eq. 

(47) reveals that J ss (nEk) and J s j (riEk) do not give a 
spin current in the weak momentum scattering regime. 
We can develop a physical understanding of this fact. 
In the absence of spin precession, skew scattering and 
side-jump scattering separate up-spins from down-spins. 
When band structure SO interactions are present, each 
spin processes about an effective magnetic field which de- 
pends on fc, thus it is not conserved. Electrons are driven 
by the external field and collide with impurities, with 
up-spins scattering predominantly in one direction and 
down-spins predominantly in the other direction. The 
spins then travel towards the edges of the sample, yet 
they are subjected to the action of the band structure 
SO effective field, which causes them to precess. Upon 
arriving at the edge the spins are completely random- 
ized. Therefore, very generally, side-jump scattering and 
skew scattering do not give rise to a spin current in 2D 
systems. 



C. Anomalous spin precession from impurities 

The last piece in the puzzle is the driving term 
J(i\(n>Ek), which needs to be studied independently for 
each model. We denote the contribution of this term to 
<7a by cr A ct . Once found, this term is added to a^ rcc to 
give oa, which yields the total SHE due to r so . 



1. Linear Rashba and Dresselhaus SO 
For linear Rashba band structure SO coupling Hr\ 
Jax{n Ek ) = E-ka-d6(k-k F ). (50) 

The spin-Hall conductivity due to this driving term is 



This term exactly cancels <7^ rec . The same holds for the 
linear Dresselhaus SO interaction Hex- 



B. Anomalous spin precession from electric field 



2. Cubic Dresselhaus SO 



Using Eq. 47, we have a term in the density matrix 



rrprec x 

S Ek\ — 2 



1 



S(e 0k - e F ). (48) 



In the weak momentum scattering limit fifcT S> 1 this 
result is independent of the form of the band structure SO 
interaction. We have given (in this subsection alone) a 
result valid beyond the weak momentum scattering limit 
so as to emphasize this apparent independence is only an 
artifact of this limit. For electron systems in this limit, 
the spin-Hall conductivity due to this term is, 



n e e\ 



(49) 



In general and c A do not cancel. We consider 

a 2DEG in a wide quantum well at high electron density 
n e , where the dominant band structure SO coupling term 
is described by the cubic Dresselhaus Hamiltonian -f/£>3. 
The scattering term J^(riEk) is given by 



t , \ meBXk 3 _ . ~ . „„ 
Jnx(riEk)= p — Ek(<r0sm26 

- er • fccos26»)5(fc - k F ). 



(52) 



This gives a trivial contribution to the spin-Hall current, 

rif.eX 



160 ' 



(53) 
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which is two orders of magnitude smaller than c^ IGC , and 
may be neglected. The remaining term due to r 80 is 
cr^ rec , and thus in the weak momentum scattering limit 



OA 



n e e\ 



(54) 



The magnitude of the SHE conductivity due to the band 
structure SO coupling (the band-structure SHE) in the 
2D cubic Drcsselhaus model has been calculated to be 
3e/87r. 69 Therefore the total SHE conductivity, including 
that due to band structure SO, is 



3e n e e\ 

h — — . 

8tt 2 



(55) 



The term due to band-structure SO is density- 
independent, whereas the anomalous spin precession 
term in the SHE is linear in n e . These are the only two 
terms, so the anomalous spin precession contribution can 
be directly measured through its density dependence. 



3. Hole systems 



It is easily seen that in 2D hole systems both o^ lec and 
tr^ ct are zero. For holes, can be found from Eqs. (14) 
and (20). Substituting this into Eq. A6, we find that the 
spin-Hall current averages to zero over directions in mo- 
mentum space. In Jn\{nEk), in ah cases studied, terms 



oc e 



±3i6 



cause the angular integral to vanish. Therefore, 



in 2D hole systems 



0. 



(56) 



There is thus no contribution to the SHE due to anoma- 
lous spin precession in 2D hole systems. 



VII. DISCUSSION 

To summarize, o~\ = in 2D hole systems, while in 
2D electron systems in the weak momentum scattering 
regime it can be written as 



n e e\ 



(57) 



The results for the total SHE due to r so are summarized 
in Table I. Interestingly, o\ can be nonzero, even though 
that is only true in one out of the several situations stud- 
ied explicitly in this work. 

We have argued previously that er^ I0C should be 
thought of not as a surviving side-jump term, but a qual- 
itatively new term, which is not present in systems with- 
out band structure SO. We demonstrate that this term 
is related to spin precession induced by both band struc- 
ture SO and r so . The electric field E gives rise to an 
additional SO effective field A& II i-direction. The band 



structure SO effective field ftk is in the plane. We ex- 
amine spin precession in the total effective magnetic field 
f2fc and Afc, redefining f2fc — > with 



fit, — fit 



(58) 



Let fife = (£l x ,0,0) and turn on E adiabatically, gener- 
ating a small Vl z <C Q x - We study the Heisenberg equa- 
tion of motion for the spin (Bloch) vector s, which reads 
ds/dt = ft x s, in a clean system. The spin is taken 
initially to be parallel to (l x . In component form 



ds x 
dt 

(lSy 

dt 

ds z 
dt 



Q Z Sy , 



(59a) 



(59b) 



(59c) 



One can take the time derivative one more time and solve 
the equations exactly, yet the physics is evident from Eq. 
(59b). Since s(t = 0) = (s x ,0,0) and s y is initially zero, 
s y should remain zero at all times. Setting ds y /dt in the 
steady state we obtain 



S 2 :r 



(60) 



The explanation is as follows: s y is initially and must 
remain 0. When E is turned on an additional component 
2 is generated, which makes s x precess and gives a small 
contribution to s y . To cancel this, s z must develop a 
small out-of-plane component, which precesses around 
Q x , and gives the exact opposite contribution to s y . The 
extra s z density has opposite signs for the two halves of 
the Fermi surface, giving rise to a net spin-Hall current. 
The argument presented here shows that r so gives rise to 
a spin-Hall current even in a clean system. We refer to 
this process as anomalous spin precession. 

This argument can be generalized to explain anoma- 
lous spin precession in a disordered system as well. This 
can be done by replacing E — ► E + VU(r), and un- 
derstanding this to represent the total local electric field. 
We thus reproduce both anomalous spin precession terms 
- the one due to the external electric field and the one 
due to the impurity potential. Both terms give an effec- 
tive magnetic field out of the plane of the quantum well, 
modifying the intrinsic SO spin precession. 

Equation (47) is valid for weak momentum scattering. 
Appendix A shows that in the strong momentum scatter- 
ing regime SEk\ diverges. Physically, this is because we 
are using Clk as our reference, and projections parallel 
and perpendicular to it become ill-defined as — > 0. 
In this limit Dyakonov-Perel spin relaxation is no longer 
active, and there is no spin relaxation at all. We demon- 
strate in Appendix A that the divergence in the strong 
momentum scattering regime is cured by the introduc- 
tion of the Elliott- Yafct spin relaxation time tey , which 
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is also related to VW. Nevertheless, in order to be con- 
sistent one would have to formulate the entire theory up 
to order A 2 , which is beyond the scope of this paper. 

In deriving a\ we have assumed for simplicity that the 
scattering potential is short ranged. We do not expect 
the results to change qualitatively for long-range impu- 
rities. Firstly, we have shown that tr^ rec is independent 
of scattering in weak momentum scattering limit and is 
traced to a mechanism unrelated to disorder. Secondly, 
for a general potential the anisotropic terms in U kk ' will 
depend strongly on the form of the potential, as will cr| ct . 
By this argument, in general we expect no cancellation 
between er A ct and (J^ Tec , even though <r| ct may have a 
different numerical value from that determined in this 
paper. Finally, past experience with the SHE shows that 
important cancellations, such as that of the SHE due to 
Rashba band structure SO coupling, tend to have a fun- 
damental origin 70 and are independent of whether the 
scattering potential is short-range or long-range. 8,9 

VIII. SUMMARY AND CONCLUSIONS 

We have determined all the contributions to the SHE 
due to the anomalous position operator f so in 2D elec- 
tron and hole systems. The SHE due to skew scattering 
and side-jump scattering vanishes in the presence of spin 
precession caused by the band structure SO coupling. 
Two additional contributions to the SHE exist due to r so , 
one of which is scattering-dependent and one of which is 



due to anomalous spin precession under the action of r so 
and the electric field. These two contributions cancel out 
in systems with band structure SO linear in k, and are 
independently zero in 2D hole systems. However, the 
contribution due to anomalous spin precession survives 
in 2D electron systems with a significant cubic Dressel- 
haus term i.e. for wide quantum wells with high electron 
densities. Anomalous spin precession can therefore be 
detected in such a system. 

A full account of the SHE in 2D systems must include 
the lengthy calculation of the electric field contribution 
to the skew scattering term, plus the band structure SO 
correction to that term. Moreover, in this work we have 
only considered heterostructures grown along the main 
crystal axes. Finally, the full answer will be known when 
the definition of the conserved spin current is taken into 
account, as has been done for the band-structure SHE. 71 
We reserve these studies for a future publication. 
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Appendix A: Decomposition of the spin density matrix into Sei^x and TEkx 



i^From the quantum Liouvillc equation, we obtain for SEkx and Teux for short-range impurities 



dSEkx , i 



dt 



— [H, Teux] = T>Ekx 



+ I [H, T Ek x] + — = (p EkX - V Ek x) - ~ [H, S Ek x] + I [H, T EkX ]. 
at n t n n 



(Ala) 
(Alb) 



On the RHS of Eq. (Alb) we substitute for - [H, T Ek x] horn Eq. (Ala). We rewrite Eqs. (Al) as 



dS 



dt 



T [H k , TE k \] — V Ek\ 



+ \ [H k ,T Ekx ] + ^ = V Ek x - ( + \ [H kl S Ek x] 



(A2a) 
(A2b) 



Defining T Ek \ = e- iHht / h f Ek \ e lHkt / h and S Ek \ = e~ %H ^/ h S Ek \ e lHkt / h , we can easily solve Eq. (A2b) 



9T Ek x T Ek x _ jHkt/h^ -iH k t/h dS Ek \ 
at t at 



TEkX = —SEk\ 



dt' e 



JHt'/H^^^-iHt'/H,, S Ek\ 



% VEk\e 



(A3a) 
(A3b) 
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where the last line was obtained by integration by parts. We can write T Ekx (without the tilde) as 



T EkX = -S^T X + J™ dt> e-^e- iHt '/ n (v EkX + ^) , 



Using S Ekx = ^ cr ■ s E k\, V Ekx = \< J ' d EkX and T E k\ = \ <J ' *sfcA, and carrying out the time integral 

S E k\ \ ^fcT 2 (d E k\T ) 



t E k\ = J~2 fe x d Ekx H 2 + — -pro— 5 + additional terms. (A5) 

\ r / 1 + fifcT 2 1 + fi 2 t 2 

The physical interpretation of the terms appearing in Eq. A5 is as follows. The first term [containing £l k x (. . .)] 
gives the full spin current when there is spin precession (D, k ^ 0). The second term (containing d Ekx r) recovers 
the spin current due to impurity SO coupling when there is no spin precession (tt k — 0). It vanishes in the weak 
momentum scattering limit D, k r 3> 1. Finally, the additional terms ensure that t Ekx averages to zero over directions 
in momentum space, but these terms give no spin current. 

Let Aij = (Sij — flki^kj), abbreviate As Ekx = AijS EkX j , and substitue Eq. (A5) into Eq.(A2a). In the steady 
state 



s EkX = d Ekx - )Ad Ekx + - ■ „ „ ' , (A6) 



1 + ft 2 r 2 



A 



i + n 2 r 2 1 + 



1 



using - [H k , S Ekx ] — — - a ■ f2 fc x s Ekx . For fi fc r > 1 we obtain simply 

A 



—j s EkX = d EkX - Ad EkX (A7a) 
t Ekx = ( — —) x [ d Ekx H — — ] + additional terms. (A7b) 

\n k j \ t 



Appendix B: Elliott- Yafet spin relaxation time 

This derivation is for a general S k . Consider the scattering term in the Born approximation Eq. (26) up to second 
order in A, and focus on its action on S k . In this term we may ignore the part of the time evolution operator oc Q k . 
This scattering term is referred to as Jey and takes the form 

JEv(S k ) = / V kk > (V k > k S k - Sfc'Vfc'fc) S(e k - ek') + h.c. (Bla) 

h J (2tt) 2 

= Tjjj^JJ Jd V ^'| 2 Sk-V kk 'Sk'Vk'k) + h.c. (Bib) 
Bearing in mind that |Vfcfc'| 2 is a scalar, and in 2D systems V kk ' oc a z , the term V kk i S k 'V k > k has two possible forms: 

f\V kk <\ 2 S k ,, foTS k >cco- z 

y kk ' ok' v k ' k — \ 

\ -\V kk >\ 2 S k ', for S k ' <xa x ,a y 



J E Y(S k ) = <E_ \ Vkk ,\^ Sk - m z S k '), (B3) 



so that 

h? J 2tt 

where m z = — 1 before o~ x , o~ y and m z = 1 before a z . If S k oc cr z the spin is out of the plane and is conserved during 
scattering, thus JEy(S k ) gives just a correction to the momentum relaxation time. The change of sign for S k oc a x , a y 
is crucial. For short-range impurities, with \V kk >\ 2 = X 2 k 4 \U\ 2 sin 2 7, 

JMS k )= — / — (5 fc -m,S fc 0(l-cos27). (B4) 

If we now write S k = S k + T k , and define (1/tey) = A 2 fc 4 /r, then JEy(S k ) simplifies to 

f /c \ S k -m z S k m z f 27r d& . . 

JEY{Sk)= s + - / — T fc , cos2 7 . (B5) 

Ztey Ztey Jq Ztt 
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Appendix C: tey cures divergence in SEk\ 

Equations (A7a) are correct as long as Sir 3> 1, otherwise s Ekx found from Eq. (A6) diverges at small fir. The 
way out of this dilemma is provided by the Elliott- Yafet spin relaxation time. Consider adding JEy(S k ) to Eqs. (Al) 



^ + - h W^\+ SEkX 2T 2 SEkX = ^ (Ola) 
dT EkX , i rTT m i , T EkX 



, -r[H,T EkX } + ^^ + J EY (T EkX ) = (V k -V h )--[H,SEkx\ + -r[H,T EhX }, (Clb) 

(71 Al T(ot II II 

where l/r tot = 1/r + 1/tey- Since A/c^ <C 1, the Elliott- Yafet spin relaxation time tey 3> t, and the term containing 
the angular integral over 9' is a very small correction to Eq. (Clb), which may be neglected. The only change to the 
above formalism is an extra term in the equation for S EkXl which is nonzero for S Ekx in plane. The spin generated 
by an electric field is in-plane, so we can focus on this component, for which m z = —1, and Eq. (A6) becomes 



i + ^fc-w T tot tey \i + n 2 k T t 2 ot j 

This cures the unphysical divergence at small fl k r. To see this, consider the simplest case, that of isotropic £l k , 



2d Ekx r tot {l + n 2 k Tl t ) - 2Ad Ekx n 2 k T? ot 

[nlrl t + (2r tot /T EY )(l + nlrl t )} ' 1 > 

Clearly s Ekx — > as Vl k — > 0. Physically, the Elliott- Yafet spin relaxation time is needed to cure this divergence 
because projections parallel and perpendicular to tt k arc ill-defined as £l k — > 0. 
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